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ABSTRACT 

We determine the radio and optical luminosity evolutions and the true distribution of the radio 
loudness parameter defined as the ratio of the radio to optical luminosity, for a set of more than 
5000 quasars combining SDSS optical and FIRST radio data. We apply the method of Efron and 
Petrosian to access the intrinsic distribution parameters, taking into account the truncations and 
correlations inherent in the data. We find that the population exhibits strong positive evolution with 
redshift in both wavebands, with somewhat greater radio evolution than optical. With the luminosity 
evolutions accounted for, we determine the density evolutions and local radio and optical luminosity 
functions. The intrinsic distribution of the radio loudness parameter R is found to be quite different 
than the observed one, and is smooth with no evidence of a bi-modality in radio loudness. The results 
we find are in general agreement with the previous analysis of Singal et al. 2011 which used POSS-I 
optical and FIRST radio data. 
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1. INTRODUCTION 

Quasars are distant active galactic nuclei (AGN) with 
emission seen across the electromagnetic spectrum, from 
radio to X-rays. The optical emission from quasars is 
thought to be dominated by the radiation from the ac- 
cretion disk around supermassive black holes, while the 
radio emission is dominated by the plasma outfiowing 
from the black hole/accretion disk systems. Because 
of this, different but complementary information can be 
gathered from both photon energy ranges regarding the 
cosmological evolution of AGN. It is therefore important 
to determine in detail the redshift evolutions of quasars 
in both radio and optical regimes. ^ 

In a previous paper (Singal eFaII[20TTI . hereafter QPl) 
we explored the luminosity evolutions and radio loud- 
ness distribution of quasars with a dataset consisting of 
the overlap of the FIRST (Faint Images of the Radio Sky 
at Twenty Centimeters) radio survey with the Automatic 
Plate Measuring Facility catalog of the Palo niar Observa- 
tory Sky Survey (POSS-I), as presented bv lWhite et al.l 
(2000). In this paper we present the results from a much 
larger dataset consisting of the overlap of FIRST with 
the Sloan Digital Sky Survey (SDSS) Data Release 7 
(jAbazajian et al.l l2009l ) quasar catalog, which contains 
a factor of ^ 10 more objects and spans redshifts from 
0.065 to 5.46. 

In the literature, the evolution of the quasar luminos- 
ity function (LF) has been described not only for opti- 
cal and radio luminosities but also for X-ray, infrared. 
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and b olometric luminosities (e.g . lUeda et aTl l2003l 
Richa rds et all 120061: iMatute et al.l 2006: H opkins et all 
2007; Groom et al.ll2009f ). The shape of the LF and its 
evolution are usually obtained from a flux limited sample 
in the a waveband fa > fm,a with fm,a denoting the flux 
limit and La = 4.7:d\{z) {1/Ka{z)) fa^ where dL is the lu- 
minosity distance and Ka{z) stands for the K-correction. 
For a pure power law emission spectrum of index Sa de- 
fined as fa cx one has Ka{z) = (1 + z)-^~^^. This 
simple form may be modifled for optical data by the pres- 
ence of emission lines, as in this work. 

In general, the determination of the LF and its evo- 
lution requires analysis of the bi-variate distribution 
^a{La^z). The flrst step of the process should be the 
determination of whether the variables of the distribu- 
tions are correlated or are statistically independent. A 
correlation between La and 2: is a consequence of luminos- 
ity evolution. In the case of quasars with the optical and 
some other band luminosity, we have at least a tri-variate 
function. One must determine not only the correlations 
between the redshift and individual luminosities (i.e. the 
two luminosity evolutions) but also the possible intrinsic 
correlation between the two luminosities, before individ- 
ual distributions can be determined. Knowledge of these 
correlations and distributions are essential for not only 
constraining robustly the cosmological evolution of ac- 
tive galaxies, but also for interpretation of related obser- 
vations, such as the extragalactic background radiation 
(e.g. Singal et al. 20_1Q|; [Hopkins et al. 2010). 

A related question is the distribution of the 'radio- 
loudness parameter', R = Lrad/Lopt^ for the quasar pop- 
ulation, defined as the ratio of the 5 GHz radio to 2500 
A optical luminosity spectral densities, and the distinc- 
tion between so-cahed 'radio loud' {R > 10; RL for 
short) and 'radio quiet' {R < 10; RQ for short) quasars. 
Weak hints of a bi-modality in the di stribution of the 
radio loundess parameter described bv iKellerman et al] 
(|1989l ) suggested that logR = 1 could be chosen as the 
radio loud/quiet demarcation value. Using this value 
for the division between RL and RQ quasars, the dif- 
ferences between the two classes have been investigated. 
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including the possibility of distinct cosmol ogical evolu- 
tion of the RL and RL populations (e.g. iMiller et al.l 
[19901 : iGoldschmidt et al.lflQOl 'Jiang et al.l 120071 ). Stih, 
the more recent analyses of different samples of objects 
reported in the literature so far gave rather inconclusive 
results on whether any bi-modality in the distribution 
of the radio loudness parameter for quasars is inherent 
in the po pulation (see'Ivezic et al."'2002'; 'Ciras uolo et al. 
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2011: Mahony et al. 2012). In QPl we found no evidence 
for a bi-modality in R in the range — 1 < logR < 4. 

In addition to the above cited works, there have been 
many papers dealing with this ratio and RL vs RQ issue, 
as well as luminosity ratios at other wavelengths, e.g. 
IR/radio, Optical/X-ray etc. However, to the best of 
our knowledge, except for QPl, none of these works have 
address the correlations between the radio and optical 
luminosities, which is necessary for such and analysis. 
Additionally, in general they have not concentrated on 
obtaining the intrinsic — as opposed to the raw observed 
— distribution (and/or evolution) of the ratio, which is 
related to the tri-variate LF ^(Lopt, i^rad, ^) b}|3 



Gr{R, z) = ^(i^opt, R ^opt, z) Lopt dLopt 
Jo 



Jo 



^ I — ^,^rad,^ I ^rad— ^2~* 



(1) 



In Appendix A of QPl we showed how, even in the 
simplest cases, the observed distributions can be very 
different from the intrinsic ones. Thus, for determina- 
tion of the true distributions the data truncations must 
be determined and the correlations between all variables 
must be properly eva l uated . 

lEfron fc PetrosianI (|1992L [1299', EP for short) devel- 
oped new methods for determination of the existence of 
correlation or independence of variables from a flux lim- 
ited and more generally truncated dataset, which were 
further expanded in QPl. Our aim in this paper is to 
take all the selection and correlation effects into account 
in determination of the true evolution of optical and ra- 
dio luminosities and their ratio and to find their distribu- 
tions, using the larger SDSS DR7 QSO x FIRST dataset. 

In ^we describe the radio and optical data used. ^ 
contains a general discussion of luminosity evolution and 
the sequence of the analysis. In g] we apply the EP 
method to achieve the luminosity-redshift evolutions and 
the correlation between the luminosities. We determine 
the density evolution in ^ and the local luminosity func- 
tions in ^ A discussion of the radio loudness distribu- 
tion is presented in ^ ^ we investigate some of the 
assumptions used, and ^ contains a discussion of the 
results. This work assumes the standard ACDM cosmol- 
ogy throughout, with Hq = 71 kms~^ Mpc~^, = 0.7 
and Qrn = 0-3. 

2. DATA 

In order to evaluate the luminosity evolution in both 
radio and optical, and to separate and compare these 

^ Equation [T] arises because by definition f GR(R,z)dR = 
f f ^(Lopt, ^rad ^) dLopt dL^^f^, and, following from the definition 
of R, dL^^^ = Lopt dR and dLopt = -(^rad/^^) dR. 




Fig. 1. — The 2500 A rest fram e absolute luminosity density 
for the SDSS DR7 quasar catalog (Schneider et aL'^MTO). Black 
points are the objects with the i < 19.1 criterion (63942 objects), 
while blue points are the objects outside of this subset (41728 ob- 
jects). It is seen that the i < 19.1 subset forms a catalog that has a 
smoother redshift distribution without a bias toward objects with 
z > 2, and with a uni form limiting fiux for every redshift (see §6 of 
ISchneider et al.l 120101 ). Also shown is the upper limiting fiux cor- 
responding to i =15.0. The 2500 A luminosity density is obtained 
from the observed i-band magnitude, converting to fiux at the inte- 
grated center band frequency, and applying the luminosity distance 
obtained from the redshift with the standard cosmology and the 
K-corrections provided by R06 which include the continuum and 
emission line effects. 
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Fig. 2. — Same as Figure [T]but for the quasars in the canonical 
SDSS X FIRST dataset used in this analysis (5445 objects). The 
black points are the RL objects while the red points are the RQ. 

effects, we require a dataset that has both radio and op- 
tical fluxes to reasonable limits across a broad range of 
redshifts that contains a significant number of both RL 
and RQ object s. The overlap of the FIRST bright q uasar 
radio survey (B ecker et aHllQQSl : iWhite et al.1ll997l) with 
the S PSS DR7 quasar catalog (DR7Q: ISchneider et al.l 
can form such a dataset. 
The SDSS DR7 quasar catalog has a limiting i band 
magnitude of 21 or fm,i = 0.015 mJy, and contains over 
105,000 objects with redshifts ranging from 0.065 to 5.46 
over 9380 deg^. We work with a more restrictive z-band 
optical magnitude limit of 19.1 {fm,i = 0.083 mJy) to 
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Fig. 3. — The 1.4 GHz rest frame absolute luminosity density for 
the quasars in the canonical SDSS x FIRST dataset used in this 
analysis (5445 objects). To obtain the 1.4 GHz luminosity density 
we use the luminosity distance obtained from the redshift and the 
standard cosmology and the standard K-correction. We assume a 
radio spectral index of 0.6. The black points are the RL objects 
while the red points are the RQ. 
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Fig. 4. — The redshift distribution of the ratio R of rest frame 
absolute luminosities at 5 GHz and 2500 A for the quasars in the 
canonical SDSS x FIRST dataset used in this analysis. The 5 
GHz luminosity is obtained from the 1.4 GHz luminosity assuming 
a radio spectral index of 0.6. The black points are the RL objects 
while the red points are the RQ. Also plotted are the observed 
fraction of objects with R > 10 (dotted curve) and R > 100 (solid 
curve) . 

form a parent optical catalog that has a smoother red- 
shift distribution with a reduced bias toward objects with 
z >2 and with a uniform limiting flux for every redshift 
(see §6 of DR7Q). The i < 19.1 catalog consists of 63,942 
objects, with a maximum redshift of 4.98, shown in Fig- 
ure [TJ 

In the SDSS DR7 catalog, objects are identified as 
quasar candidates if either they have the requisite opti- 
cal colors, or if they have a radio match within 2''. Only 
1% of the quasars in the catalog were selected from the 
candidate list for spectroscopic followup based on the 
later criterion. We also note the presence of an upper 
limit i band magnitude of 15.0 for inclusion in the cat- 
alog. However this criterion is not completely rigorous. 



as mentioned in DR7Q. Further properties of the DR7 
quasar catalog, such as black hole masses ob tained from 
l ine em ission measurements, are presented in lShen et al.l 
(|20TT1 ). 

The FIRST survey, carried out with the Very Large Ar- 
ray (VLA) in B configuration between 1995 and 2004, has 
a limiting peak pixel 1.4 GHz fiux of 1 mJy, and contains 
816,000 sources over 9500 deg^. Different criteria can be 
used to determine radio and optical matche s to form a 
joint SDSS DR7 QSO x FIRST catalog. 'Jiang et_alj 
(|2007[ ). hereafter J07, in their analysis of SDSS quasars, 
used a matching radius of 5'' for single radio matches 
and 30'' for multiple matches to an optical source (al- 
though they included quasars with no radio detection in 
their joint sample while we do not). We construct a fidu- 
cial catalog using the radius matching criteria, but, as 
mentioned, unlike J07 we do not includ e qu asars with no 
radio detection (see the discussion in § 19. ip . 

This fiducial catalog contains 5677 objects ranging in 
redshift from 0.064 to 4.82 and with \ogR ranging from 
-0.47 to 4.44. 4327 of the objects in the fiducial catalog 
are single matches (which J07 calls 'Fanar off- Riley type 
I sources', FRls, perhaps incorrectly; see the discussion 
in §9.ip and 1350 are multiple matches (which J07 calls 
'Fanar off- Riley type II objects', or FR2s for short). This 
sample seems to be skewed toward RL as opposed to RQ 
quasars, with 4543 (80%) having R > 10. It should be 
noted that J07 used the SDSS Data Release 3, so our 
realization with the same matching criteria, even with 
the additional i < 19.1 magnitude cut and no objects 
without radio detections, results in more objects (J07 
have 2566 objects). 

We also construct catalog with a universal 5'' matching 
criterion which we feel is more appropriate, as discussed 
in ^9.11 In this catalog there are 5445 objects with only 
25 multiple radio matches to an optical source. We will 
refer to this as the canonical sample, and base our analy- 
sis primarily on it, while comparing to results we obtain 
with the fiducial sample, which as discussed below are 
quite similar. In the canonical sample 4466 (82%) have 
R > 10. The major changes from the fiducial catalog to 
the canonical one are that some optical objects with mul- 
tiple radio matches in the fiducial sample are reduced to 
having only single matches in the canonical sample, with 
a corresponding reduction in radio fiux, while a small 
number of objects drop out entirely, having had multi- 
ple matches within 30'' but none within 5". The fiducial 
sample and canonical sample are nearly identical in their 
redshift distribution, and the average radio luminosity 
differs only slightly, being 10% lower in the canonical 
sample. 

As shown in [Richards et al.l (j2006 f). hereafter R06, e.g. 
Figure 8 of that work, the SDSS optical quasar sam- 
ple contains significant biases in the redshift distribution 
due to emission line effects and the differing fiux limits 
at z > 2. We have reduced the later effect by using only 
the universally fiux-limited i < 19.1 sample, and the for- 
mer by adopting the full ^opt = 0.5 power law continuum 
plus emission line K corrections presented in R06 and dis- 
cussed in §5 and Table 4 of that work. The methods of 
this work can then account for any bias resulting from 
emission line effects, as long as they are included in the 
conversions from luminosity to fiux, i.e. in the K correc- 
tions. For the radio data we assume a power law with 



4 



Singal et al. 



^rad = 0.6. Figures [2] and [3] show the radio and optical 
luminosities versus redshifts of the quasars in the canon- 
ical constructed SDSS x FIRST catalog. Figure |4] shows 
the radio loudness parameter R versus redshift for the 
canonical dataset. 

3. GENERAL REMARKS ON LUMINOSITY FUNCTIONS 
AND EVOLUTIONS 

3.1. Luminosity and density evolution 

The LF gives the number of objects per unit comoving 
volume V per unit source luminosity, so that the number 
density is dN/dV = J dL^^ ^^{La^ z) and the total num- 
ber is A/" = J / dz {dV / dz)^ Q,{La^ z). To examine 
luminosity evolution, without loss of generality, we can 
write a LF in some waveband a as 



^a{La, Z) = p{z) i)a{Lal9ciz), rf^lg^z), 



(2) 



where gdz) and describe the luminosity evolution 
and comoving density evolution with redshift respec- 
tively and stands for parameters that describe the 
shape (e.g. power law indices and break values) of the 
a band LF. In what follows we assume a non-evolving 
shape for the LF (i.e. rf^ — const, independent of L and 
z\ which is a good approximation for determining the 
global evolutions. Once the luminosity evolution gdz) 
is determined using the EP method we can obtain the 
mono-variate distributions of the independent variables 
L'^ = La/^a(^) and z, namely the density evolution p{z) 
and "local" LF i/jg,. The total number of observed objects 
is then 



tot 



lo i Lr 



dV IpaiLa/gaiz)) 

dLa p[z) , 

dz ga{z) 



(3) 

We consider this form of the LF for luminosities in 
different bands, allowing for separate (optical and radio) 
luminosity evolutions. 

In the previous work (QPl) with the White et al. 
dataset, we assumed a simple power law for the lumi- 
nosity evolution 



g,{z) = {l + zY 



(4) 



However with the SDSS dataset which contains objects 
past z=4 and many more objects past z = 3, we need to 
use a more complicated parameterization: 



5'a(^) = 



(1 + ^)^- 



(5) 



With either definition of ga{z) which gives ^a(O) = 1, the 
luminosities LJ^ refer to the de-evolved values at z = 0, 
hence the name "local" . 

We have verified the form of equation [5] as a good fit 
to the evolution by considering the optical evolution in 
the large parent optical only sample in narrow redshift 
bins with the appropriate simple form of equation |4] and 
then verifying the form of the cumulative evolution over a 
larger range. We have determined the appropriate value 
for Zcr as well as confirmed that the same exponent is ap- 
propriate in the numerator and denominator, using the 
large set of optical only data and the methods of §4.11 
With the exponents fixed to be the same, the optimal 



value of Zq^ is determined to be 3.7±0.3, by considering 
the Icr range of the best fit evolution while letting that 
numerical factor vary. We would expect that if the trun- 
cations and correlations have been properly accounted 
for by the methods of this work, ^opt(^) as determined 
from the simultaneous radio and optical dataset should 
match that as determined from the parent optical only 
dataset, which is shown to be the case (see 

We discuss the determination of the evolution factors 
ga{z) with the EP method, which in this parameteriza- 
tion becomes a determination of /c^, in 21 The den- 
sity evolution function p{z) is determined by the method 
shown ^ Once these are determined we construct the 
local (de-evolved) LF ip^^L'^)^ shown in ^ 

3.2. Joint Luminosity Functions 

In general, determination of the evolution of the LF 
of extragalactic sources for any wavelength band except 
optical involves a tri-variate distribution because redshift 
determination requires optical observations which intro- 
duces additional observational selection bias and data 
truncation. Thus, unless redshifts are known for all 
sources in a radiol survey, we need to determine the com- 
bined LF ^(Lopt, ^rad, ^) from a tri-variate distribution 
of z and the fiuxes in the optical and radio bands. If 
the optical and radio luminosities were statistically in- 
dependent variables, then this LF would be separable 

^(Lopt,i^rad,^) = ^opt(^opt,^) X ^rad(^rad,^) and WC 

would be dealing with two bi-variate distributions. 

However, the luminosities in the different wavebands 
are correlated. The degree and form of the correlation 
must be determined from the data, and as described be- 
low, the EP method allows us to determine whether any 
pair of variables are independent or correlated. Once it 
is de termined that the luminosities are correlated (see 
§4.2p . the question must be asked whether this luminos- 
ity correlation is intrinsic to the population, or induced 
in the data by fiux limits and/or similar luminosity evo- 
lutions with redshift. Determination of which is the case 
is quite intricate as discussed in Appendix B of QPl, and 
has not been explored much in the literature. This will 
be the subject of a future work. Here we will simply 
consider both possibilities. 

At one extreme, if the luminosity correlation is intrin- 
sic and not induced, one should seek a coordinate trans- 
formation to define a new pair of variables which are 
independent. This requires a parametric form for the 
transformation. One can define a new luminosity which 
is a combination of the two, which we call a "correlation 
reduced radio luminosity" Lcrr = ^rad/^(^opt/^fid)5 
where the function F describes the correlation between 
Lrad and Lopt and Lfid is a fiducial luminosity taken here 
to be 10^^ ergsec"^ Hz~^ For the correlation function 
we will assume a simple power law 



Lc 



(iopt/ifld)" 



(6) 



where a is a bulk power law correlation index to be de- 
termined from the data. This is essentially a coordinate 

^ This is a convenient choice for as it is lower than the 

lowest 2500 A luminosity considered in our sample, but results do 
not depend on the particular choice of numerical value. 
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rotation in the log- log luminosity space. As shown in §H 
below, EP also prescribe a method to determine a best 
fit value for the index ol which orthogonalizes the new 
luminosities. Given the correlation function we can then 
transform the data (and its truncation) into the new in- 
dependent pair of luminosities (Lopt and Lcrr)- 

At the other extreme, if the correlation between the 
luminosities is entirely induced, then the LPs are sepa- 
rable as described above and the analysis can proceed 
from there. As mentioned, we will consider both pos- 
sibilities here. It turns out that the results obtained in 
both cases are very similar, implying that the luminosity- 
redshift correlations in the analyzed sample are much 
stronger than any intrinsic luminosity-luminosity corre- 
lations (if present). Purther mechanics of the procedure 
for intrinsically correlated luminosities is discussed in the 
Appendix. 

In this work we are also interested in the distribu- 
tion of radio loudness parameter i?, specifically its de- 
evolved value B! = L'^^JL'^^^. In the case of uncorre- 
lated radio and optical luminosities or an induced, non- 
intrinsic radio-optical luminosity correlation, the local 
distribution of B! ^ denoted hereafter as G{B!\ can be 
constructed from the local optical and radio luminosity 
functions as in equation [TJ with the redshift evolution 
function 



In the case of an intrinsic radio-optical luminosity corre- 
lation, Gr/ and the evolution can be constructed by 
equations [3] and H] given in the Appendix. 

4. DETERMINATION OF BEST FIT CORRELATIONS 

We now describe results obtained from the use of the 
procedures described in ^ on the data described in ^ 
Here we first give a brief summary of the algebra involved 
in the EP method. We generally follow the procedures 
described in more detail in QPl. This method uses the 
Kendall tau test to determine the best-fit values of pa- 
rameters describing the correlation functions using the 
test statistic 

x = 5i^=£d (S, 

to test the independence of two variables in a dataset, 
say (xj^i/j) for j = 1, . . . , n. Here 1Zj is the dependent 
variable {y) rank of the data point j in a set associated 
with it. Por a untruncated data (i.e. data truncated 
parallel to the axes) the set assocated with point j in- 
cludes all of the points with a lower independent variable 
value {xk < Xj). If the data is truncated one must form 
the associated set consisting only of those points of lower 
independent varialbe (x) value that would have been ob- 
served if they were at the x value of point j given the 
truncation. As an example, if we have one sided trun- 
cations as in Pigures [2] and [3] — ignoring the upper op- 
tical fiux limit for the moment, then the associated set 
Aj = {k : yk > yj, y^ < yj}, where y~ is the limiting 
y value of data point j (see EP for a full discussion of 
this method). 




3.0 3.2 3.4 3.6 3.8 4.0 

optical k 

Fig. 5. — The value of the r statistic as given by equation [8] as 
a function of /copt for the form of the optical luminosity evolution 
given by equation [S] for the 63,942 quasars in the i < 19.1 optical 
only dataset. The 1 a range for the best fit value of a is where 
|r|<l. 

If (xj^yj) were independent then the rank IZj should 
be distributed uniformly between and 1 with the ex- 
pectation value and variance Ej = (1/2) (n + 1) and 
Vj = (1/12) (n^ + 1)7 respectively, where n is the number 
of objects in object j's associated set. Independence is 
rejected at the m a level if | r | > m. To find the best 
fit correlation the y data are then adjusted by defining 
Vj — Vj/^i^j) rank test is repeated, with differ- 

ent values of parameters of the function F. 

4.1. Optical- Only Dataset Luminosity- Redshift 
Correlation 

With the optical only dataset (here we use the i < 19.1 
set described in determination of the luminosity- 
redshift correlation function ^opt(^) by equation [5] reduces 
to determination of the value of the index /Copt- As ev- 
ident from Pigure [T] the Lopt — z data are heavily trun- 
cated due to the fiux limits of SDSS. The associated set 
for object j includes only those objects that are suffi- 
ciently luminous to exceed the opitlcal ffux minimum for 
inclusion in the survey if they were located at the redshift 
of the object in question, ie Lk > Lmin{zj)- As the val- 
ues of kopt are adjusted and Lopt is scaled by gopt{k,z), 
the luminosity cutoff limits for a given redshift are also 
scaled by gopt{k,z). 

Pigure [S] shows the value of the test statistic r for 
a range of values of /Copt given the parameterization of 
equation m It is seen that the best fit value of /Copt is 3.3 
with a Icr range of 0.1. This indicates that quasars have 
strong optical luminosity evolution, i.e. the LP shifts to 
higher luminosities at larger redshift s. 

We have performed an optimization procedure to de- 
termine the best-fit value of Zcr, the numerical constant 
in the denominator of equation [5l as well as whether the 
same exponent is appropriate for the numerator and de- 
nominator. This was done by varying these parameters 
and determining the best-fit range of values in the same 
way as described above. We find Zqy — 3.7±0.3. 

In this analysis we have ignored the upper optical fiux 
limit of SDSS quasars corresponding roughly to a z band 
magnitude of 15.0 as discussed in ^ The reason for this 
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Fig. 6. — The value of the r statistic as given by equation [8] as a 
function of ol for the relation Lrad oc (Lopt)"", where Lopt and Lrad 
are the optical and radio luminosities, respectively, for the quasars 
in the canonical combined dataset. The Icr range for the best fit 
value of OL is where | r | < 1. It is seen that the observed optical 
and radio luminosities are strongly positively correlated, with a 
linear relation, although t his may not be the intrinsic correlation, 
as discussed in ^4.21 ^3.21 and the Appendix. 

is that data truncations are only consequential in this 
analysis if the truncation is actually denying the sample 
data points that would otherwise be there. As can be 
seen in Figure [TJ there are very few objects approaching 
this upper truncation limit, indicating that the trunca- 
tion is not consequential, i.e. it does not appreciably 
alter the sample from the underlying population. 

4.2. Radio- Optical Luminosity Correlation 

Turning to the combined radio and optical dataset, we 
will determine the correlation between the radio and op- 
tical luminosities. The radio and optical luminosities are 
obtained from radio and optical fluxes from a two-flux 
limited sample so that the data points in the two di- 
mensional flux space are truncated parallel to the axes 
which we consider to be untruncated. Since the two lumi- 
nosities have essentially the same relationship with their 
respective fluxes, except for a minor difference in the K- 
correction terms, we can consider the luminosity data 
to also be untruncated in the Lopt-i^rad space. In that 
case the determination of the associated set is trivial and 
one is dealing with the standard Kendall tau test. As- 
suming the correlation function between the luminosities 
F{x) = we calculate the test statistic r as a function 
of a. Figure [6] shows the absolute value of the test statis- 
tic r vs a, from which we get the best fit value of a = 1.2 
with one a range ±0.1. The result similar whether the 
canonical or fiducial combined dataset is used. As ex- 
pected a is near unity, and this result is compatible with 
that obtained in QPl {a = 1.3 ± 0.2). 

As discussed in §3.2| this correlation may be inherent 
in the population or may be an artifact of the flux limits 
and wide range of redshifts. The general question of 
determining whether an observed luminosity correlation 
is intrinsic or induced will be explored in a future work. 
For this analysis going forward, we will consider both 
possibilities. We see that it does not make a significant 
difference for the major conclusions of this work. 



4.3. Joint Dataset Luminosity- Redshift Correlations 

The basic method for determining simultaneously the 
best fit kopt and /Crad, given the evolution forms in equa- 
tionO is the same as described in §4.1l but in this case the 
procedure is more complicated because we now are deal- 
ing with a three dimensional distribution (Lrad, ^opt, ^) 
and two correlation functions (^rad(^) and ^opt(^))- 

Since we have two criteria for truncation, the associ- 
ated set for each object includes only those objects that 
are sufficiently luminous in both bands to exceed both 
flux minima for inclusion in the survey if they were lo- 
cated at the redshift of the object in question. Conse- 
quently, we have a two dimensional minimization prob- 
lem, because both the optical and radio evolution factors, 
9opt{z) and ^rad(^), come into play, as the luminosity cut- 
off limits for a given redshift are also adjusted by powers 

of /Copt and /Crad- 

We form a test statistic Tcomb = \J ^opt + ^rad where 
Topt and Trad are those evaluated considering the objects' 
optical and correlation reduced radio luminosities, re- 
spectively. The favored values of /Copt and krad are those 
that simultaneously give the lowest Tcomb and, again, we 
take the la limits as those in which Tcomb < 1- For visu- 
alization. Figure [3 shows a surface plot of Tcomb- Figure 
[8] shows the best fit values of kopt and kr^d taking the 1 
and 2 a contours, for the canonical radio-optical dataset. 
We have verified this method with a simulated Monte 
Carlo dataset as discussed in QPl. 

We see that positive evolution in both radio and opti- 
cal wavebands is favored. The minimum value of Tcomb 
favors an optical evolution of kopt = 3.5 and a radio evo- 
lution of krad = 5.5, with a very small uncertainty at 
the Icr level. The fact that /Copt as determined here from 
the combined dataset is equal to that determined from 
the much larger optical only dataset (see §4.11 and Fig- 
ure [5]) indicates that the truncations have been properly 
handled. 

In the above analysis we have assumed sharp trunca- 
tion boundaries and that the data is complete above the 
boundaries. As discussed in ^this may not be the case 
for the FIRST radio data. If we restrict the sample to 
only data with /^ad > 2 mJy the favored optical evo- 
lution range for the canonical dataset lowers slightly to 
^opt = 3.0 ± 0.5 and the favored radio evolution range 
lowers slightly to k^ad = 5.0 ± 0.25. These overlap with 
the results considering the entire combined canonical 
sample at the 1 a level. 

These results indicate that quasars have undergone a 
significantly greater radio evolution relative to optical 
evolution, in general agreement with the result in QPl, 
although the best-fit differential evolution was even more 
dramatic there, given the form of g{z) used. 

The results for the best fit /Copt and k^ad are identical if 
the different fiducial radio-optical dataset matching ra- 
dius criteria are used, indicating that the addition of ex- 
tra radio flux and a few additional sources do not matter 
for this determination. If we consider that the radio- 
optical luminosity correlation is inherent in the data, 
then the orthogonal luminosities are Lopt and Lcrr (see 
§3.2p and we can determine the best fit evolutions gcv^z) 
and ^opt(^)- These results favor /copt = 3.5 and kcw = 2. 
In this case the best fit radio evolution can be recovered 
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Fig. 7. — Surface plot of the value of Tcomb for the canonical 
combined radio-optical dataset as a whole showing the location of 
the minimum region where the favored values of kopt and kcrr lie. 
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Fig. 8. — The la and 2cr contours for the simultaneous values of 
fcopt and fcrad for the canonical combined dataset, for the forms of 
the luminosity evolutions given by equation [5] 

by equation[2]in the Appendix and it is /crad = 5.5, in per- 
fect agreement with the results obtained from considering 
kept and k^ad as orthogonal, indicating the robustness of 
the result. 

5. DENSITY EVOLUTION 

Next we determine the density evolution p{z). One can 
define the cumulative density function 



''^'^ = 1 f ^^^^^ 



(9) 



which , fohowing i Petrosiaiil ()1992h based on iLvnden-Belll 

can be calculated by 



<T(.)=n(i 



m(j) 



(10) 



in object j^s associated set. In this case, the associated 
set is again those objects with sufficient optical and ra- 
dio luminosity that they would be seen if they were at 
object j's redshift. The use of only the associated set for 
each object removes the biases introduced by the data 
truncation. Then the density evolution p{z) is 



da{z) 
dz 



1 



dV/dz 



(11) 



However, to determine the density evolution, the pre- 
viously determined luminosity evolution must be taken 
out. Thus, the objects' optical and radio luminosities, 
as well as the optical and radio luminosity limits for in- 
clusion in the associated set for given redshifts in the 
calculation of a by equation [TOl are scaled by taking out 
factors of gopt{z) and ^rad(^), determined as above. Fig- 
ure [9] and [TDl show the cumulative and differential density 
evolutions, respect ivelyH The number density of quasars 
seems to peak at between redshifts 1 and 1.5, earlier 
than generally th ought for the most luminous quasars 
(e.g. iShaver et al.|[l996i) . and earlier than that found in 
R06, but more similar to the peak f ound for less lumi- 
nous quasars by Hopkins et al. (2007), and in agreement 
with Maloney & Petrosian (1999). This is slightly earlier 
than we found in QPl. 

The normalization of p{z) is determined by equation 
[31 with the customary choice of ip{L') dL' = 1. As 

stated in R06, the main sources of bias in the redshift 
distribution of the SDSS quasar sample are 1) the differ- 
ing magnitude limits for z > 2, 2) the effects of emission 
lines on i band ffux at different redshifts, and, at a some- 
what less important level, 3) the inclusion of extended 
sources at the lowest redshifts {z ~ 0.3). As discussed in 
^ we have, following R06, dealt with issue 1 by restrict- 
ing the sample to a universal i < 19.1 magnitude limit, 
and issue 2 by adopting the R06 K corrections which in- 
clude the effect of emission lines as well as the continuum 
spectrum. Since we do not see any hint of a peak or de- 
viation in p{z) or a{z) at z ~ 0.3, we do not deem issue 3 
to significantly effect our conclusions. The density evolu- 
tion p{z) can be fit with a broken polynomial in with 
coefficients shown in Table 1. 

Knowing both the luminosity evolutions gg[z) and the 
density evolution p{z), one can form the luminosity den- 
sity functions £a{z) which are the total rate of produc- 
tion of energy of quasars as a function of redshift, by 



I 



dL^ La p{z) ^dz =<L'^> ga{z) p{z) 



dV 

dz 
(12) 

We show both the 2500 A and optical and 1.4 GHz radio 
luminosity densitiy functions in Figure [TTl As expected, 
because of the greater luminosity evolution in the radio 
band, the radio luminosity density peaks earlier in red- 
shift than the optical luminosity density. 

6. LOCAL LUMINOSITY FUNCTIONS 
6.1. General Considerations 



where j runs over all objects with a redshift lower than 
or equal to and m{j) is the number of objects with 
a redshift lower than the redshift of object j which are 



^ We note a notional difference with QPl whe re w e plotted 
p(z) = da/dz rather than p as defined in equation 1111 here with 
a factor of dV/dz taken out. 
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Fig. 9. — The cumulative density function cr(z) vs. redshift for 
the quasars in the canonical dataset. The normalization of cr(z) is 
determined as described in ^polynomial fit to p{z) is given there. 
A polynomial fit to (t{z) is used to determine f(^z) by equation 1111 



0.0006 



0.0005 



0.0004 



0.0003 



0.0002 



0.0001 



0.0000 





'1 1 1 1 ' ' 
















: ^ 












'- ^ 






: ^ 






'- yii 


M 










: )K 






















: ^ 


























, *l 





2 

redshift 



Fig. 10. — The density evolution p(z) vs. redshift for the for 
the quasars in the canonical dataset. p(z) is defined such that 
ct{z) = Jq°° p(z) dV/dz dz. The normalization of p(z) is determined 
as described in JS] and polynomial fits of p(z) to z are given there. 
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Fig. 11. — The optical (solid curve) and radio (dashed curve) lu- 
minosity densities £opt{z) and £rad(z)^ discussed in ^and equa- 
tion^] Note the different scales. 



TABLE 1 
Coefficients for polynomial fit^ 

TO DENSITY EVOLUTION p(z) VS. Z 





z < 1.3 


z > 1.3 


c 


-0.000581965 


0.00487286 


ai 


+0.00752870 


-0.0145420 


a2 


-0.0233918 


+0.0196557 


as 


+0.0355617 


-0.0138312 


a4 


-0.0257837 


+0.00544187 


as 


+0.00714230 


-0.0012073 


ae 





+0.000140803 


ar 





-6.68658e-006 



^ Polynomial fits are of the form p{z) = 
ci + ai 2; + 02 2^ + as + 04 2^ + 
05 -\- ae z^ + aj z'^ . 

In a parallel procedure we can use the local (redshift 
evolution taken out, or 'de-evolved') luminosity distri- 
butions (and de-evolved luminosity thresholds) to deter- 
mine the 'local' LFs ipa{Lj)^ where again the subscript a 
denotes the waveband, and the prime indicates that the 
luminosity evolution has been taken out. We first obtain 
the cumulative LF 



(13) 



which, following iPetrosianI (|1992[ ), ^^(^a)' t)e calcu- 
lated by 



$„(L;)=n(i 



n{k) 



(14) 



where k runs over all objects with a luminosity greater 
than or equal to L^, and n(k) is the number of objects 
with a luminosity higher than the luminosity of object k 
which are in object /c's associated set, determined in the 
same manner as above. The luminosity function ipJyL'^) 
is 



dL' 



(15) 



In 0we determined the luminosity evolution for the 
optical luminosity Lopt and the radio luminosity L^ad- 
We can form the local optical ?/^opt(^opt) radio 
'0rad(^rad) Straightforwardly, by taking the evolu- 
tions out. As before, the objects' luminosities, as well 
as the luminosity limits for inclusion in the associated 
set for given redshifts, are scaled by taking out factors of 
^rad(^) and ^opt(^), with k^^d and kopt determined in gl 
We use the notation L ^ L' = L/g{z). For the local lu- 
minosity functions, we use the customary normalization 
J^/ iIj{L') dV = 1. This normalization may be biased 

by around 8% due to quasar variability as discussed in 



6.2. Local optical luminosity function 
Figures [T2l and [T3l show the local cumulative <l>opt(^opt) 
and differential ?/^opt(^opt) local optical LFs of the quasars 
in the canonical combined dataset, and that of the optical 
only dataset. 

The optical LF shows possible evidence of a break at 
~ 10^0 erg see" ^ Hz~^, which was present already in data 
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analyzed in iPetrosianI (|1973[ ) and was also seen in QPl. 
With the optical only dataset, fitting a broken power law 
yields values — 2.8 ± 0.2 and — 4.1 ± 0.4 below and above 
the break, respectively. With the canonical combined 
radio-optical dataset, fitting a broken power law yields 
values —2.8 ± 0.3 and —3.8 ± 0.5. If we allow for the 
possibility of additional uncertainty resulting from the 
consideration of possible radio incompleteness at faint 
fluxes (see discussion in the range on the power law 
below the break changes to — 2.1 ± 0.1, and above the 
break the results are not affected. 

As the optical LF has been studied extensively in vari- 
ous AGN surveys, we can compare the slope of ?/^opt(^opt) 
obtained here to values reported in the literature. In 
QPl, we found power law values of —2.0 ± 0.2 and 
—3.2 ± 0.2 below and above the break, respectively. 
iBovle et al.l (|2000D, using the 2dF optical dataset (but 
with no radio overlap criteria) use a customary broken 
power law form for the LF, with values ranging from 

— 1.39 to —3.95 for different realizations, showing rea- 
sonable agreements 

6.3. Local radio luminosity function 

Figure [Ml shows the local radio LF ipradiL'rad) - -"-^ 
seen that the local radio LF contains a possible break 
around 2 x 10^^ ergsec"^ Hz~^, with a power law slope 
of —1.5 ± 0.1 below the break and —2.6 ±0.1 above the 
break. The range for the power law above the break is 
increased slightly to — 2.5 ± 0.3 if the effects of possible 
radio incompleteness are included, as in ^ while below 
the break it is unchanged, indicating this is a small ef- 
fect. It is interesting to note that the value of the radio 
loudness parameter R corresponding to the break in the 
radio and optical luminosity functions is very close to the 
critical value ~ 10 widely discussed in the context of 
the RL/RQ dichotomy. 

We can also construct the local radio LF from 
'0opt(^opt) V^crX^crr) with cquatiou [T] of the Ap- 
pendix, under the assumption that the radio and optical 
luminosity correlation is intrinsic. This is shown by the 
small black points in Figure [HI It is seen that '0rad(^rad) 
constructed in this way is nearly identical to that deter- 
mined by considering the radio and optical luminosity 
correlation is induced, indicating that the radio luminos- 
ity function result is robust to this consideration. 

In QPl for the local radio luminosity function we found 
a break at 10^^ ergsec"^ Hz~^, with a power law slope of 

— 1.7 ±0.1 below the break and —2.4 ±0.1 above it, in ex- 
cellent agreement with the results here. The slope above 
the br eak s een here is s imilar to earlier results of Schmidt! 
(jl972[ ) and IPetrosianI (jl973[ ) which probed only high lu- 
minositie s. A more complet e comparison can be made 
with Ma uch fc Sadlerl (j2007[ ). who form radio LFs of lo- 
cal sources in the Second Incremental Data Release of 
the 6 degree Field Galaxy Survey (6dFGS) radio catalog. 
For the sources they identify as AGN, they find a break 
at 3.1 X 10^^ ergsec-^ Hz-\ with slopes of -2.27 ± 0.18 
and —1.49 ± 0.04 above and below the break, in good 
agreement with our results. 

^ It should be noted that they parameterize evolution differently 
and work in absolute magnitudes rather than luminosities, how- 
ever the slopes of their fits to the LF as they parameterize it are 
applicable, as can be seen in their section 3.2.2. 
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Fig. 12. — The cumulative local optical luminosity function 
<l>opt(^opt) quasars in the canonical combined dataset 

(stars) and for the parent optical only dataset (squares). A piece- 
wise quadratic fit to ^L'^^^) is used to determine '0opt(-^opt) 
equation 1151 The normalization of the optical luminosity function 
may be biased by around 8% (see ^) 
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Fig. 13. — The local optical luminosity function '0opt(i^opt)- 
squares show results for the parent optical only dataset, while the 
larger black stars, smaller red stars, and smaller blue stars show the 
results for the entire canonical combined radio-optical dataset, and 
that dataset with imposed radio fiux limits of 2 mJy and 4 mjy, 
respectively. The normalization of the local luminosity functions 
is described in JH] 

In Figure [TH we also show the radio luminosity func- 
tions for the same range of L'^^d^ determined by 
iKimball et al.l (|2011[ ) for two samples, one for a sample of 
SDSS optical and NRAO Very Large Array Sky Survey 
(NVSS) radio data, and another for deep radio obser- 
vations with the Extended Very Large Array in which 
almost every SDSS quasar in a narrow redshift range in 
a field was detected in radio. The results presented there 
are the luminosity function only for sources in the red- 
shift range 0.2 < z < 0.3, which we have scaled by means 
of equation [5] and converted from 6 GHz to 1.4 GHz with 
a radio spectral index of 0.6 to obtain the local luminos- 
ity function at 1.4 GHz for direct comparison with our 
results. It is seen that our luminosity functions agree on 
the faint end of the range considered here but potentially 
disagree only on the very bright end. 
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Fig. 14. — The local radio luminosity function '0rad(^rad) 
quasars in the canonical combined radio-optical dataset. The larger 
black stars show the results from considering the entire canonical 
combined dataset and the radio-optical luminosity correlation to 
be entirely induced, while the small stars show the results from the 
entire canonical combined dataset considering the correlation to 
be entirely intrinsic, calculated with equation ^ in the Appendix. 
The smaller red and blue stars show the results from considering 
an induced correlation with the canonical combined dataset and 
imposed radio flux limits of 2 mJy and 4 mJy respectively. It is 
seen that the local radio luminosity function is not sensitive to the 
assumption of whether t he ra dio-optical luminosity correlation is 
intrinsic or induced (see ^3.21 for a discussion). The normalization 
of the local luminosity functions is described in JH] We also over- 
plot the local luminosity functions as determined by Kimball et alj 
(2011) with an SDSS x NVSS sample (black crosses) and by a deep 
EVLA sample (dark blue crosses), which have been properly scaled 
for overplotting here as discussed in ^6.31 

7. DISTRIBUTION OF RADIO LOUDNESS RATIOS 

As stated in the introduction, naively one may expect 
that because the ratio R is independent of cosmological 
model and nearly independent of redshift, the raw ob- 
served distribution of R would provide a good represen- 
tation of the true distribution of this ratio. In Appendix 
A of QPl we show that this assumption would likely be 
wrong. In Figure [15] we show this raw distribution of 
R by the triangles, arrived at by using the raw values 
of R from the data and forming a distribution in the 
manner of equations [T4l and [T5l with no data truncations 
accounted for. As discussed in ^ we can reconstruct 
the local distribution of GRf{R')^ as in equation [TJ which 
provides for a more proper accounting of the biases and 
truncations and is an estimate of the true, intrinsic dis- 
tribution. The results of this calculation are also shown 
in Figure [15] by the filled circles. We include a range of 
error for the extremal possibilities that the correlation 
between the radio and optical luminosities is entirely in- 
duced or entirely intrinsic (open circles). In the later 
case, Gr'{R') is given by equation [3] of the Appendix. 

The distribution Gr'{R') calculated from the data tak- 
ing into account the evolutions and truncations is clearly 
different than the raw distribution, and shows no evi- 
dence of bi-modality in the range of R considered. The 
raw distribution is seen to be weighted toward much 
higher values of R. Using the naive, raw distribution 
instead of the reconstructed intrinsic one would result in 
estimates of quasar properties that would be significantly 
and systematically biased. 

We also know that the best fit redshift evolution of the 



ratio, given equation [4] The change in the distribut ion of 
R with increasing redshift is also shown in Figure I15F^ 
Another way to look at this is that we have found that 
the radio luminosity evolves at a different rate than the 
optical luminosity, with the consequence that their ratio 
is a function of redshift. The radio loudness of the pop- 
ulation increases by a factor of 3 by redshift 1, a factor 
of 8 by redshift 3, and a factor of 11 by redshift 5. This 
is a less dramatic evolution of the ratio than we found 
in QPl, where we used a simpler parameterization (see 
^3.ip , although still quite significant. The general trend 
toward increasing R with redshift can be validated in 
a simple way by examining the median value of R vs. 
redshift in bins of redshift which indeed shows a steady 
increase. 

This differential evolution is in disagreement with the 
result presented by JOT who show a decrease in fraction of 
RL sources with increasing redshift, which could be the 
case if the radio luminosity were to evolve more slowly 
than the optical luminosity. They however do not deter- 
mine indiv idual evolutions or LFs. On the other hand, 
iMiller et al.l (^990) have noted that the fraction of RL 
quasars may increase with redshift, which they attribute 
to a difference in the evoluti o ns of the two populations 
(RL and RQ). [D onoso et al. l ('2009') compute radio and 
optical L Fs at different re dshifts and reach the same con- 
clusion. iCirasuolo et al.l (|2Q06i) also find that the radio 
loud fraction may modestly increase at high redsh i ft. Al- 
though not directly comparable, L aFranca et al.l (|2010l ) 
show a similar evolution for R^^ the ratio of radio to X- 
ray luminosity, as we show here for R. As discussed in 
§9.H we believe that our differences with JOT stem from 
our inclusion of only radio-detected quasars in the joint 
radio-optical sample. 

We note that our results favor one population in the 
range of R considered here, in the sense that the dis- 
tribution of G{R)^ recovered from considering the data 
truncations inherent in the survey and correlations be- 
tween the luminosities, is continuous and without any 
bi-modality. This is in agreement with the conclusion 
we reached in QPl. One may hypothesize about a bi- 
modality that is not centered on the commonly assumed 
value of R= 10 or above that value, but at a lower value 
of R below those values probed in this analysis. Our re- 
sults disfavor that as well. As seen in Figures [13] and 
[EJ the lack of a significant change in the shapes of the 
computed local optical or radio LFs when the radio fiux 
limits for inclusion in the analysis are changed, or the 
parent optical only dataset is considered, argues against 
the presence of an additional population of quasars which 
would become more prominent below the lowest values of 
R considered here. The results of Kimball et al. (20iT|), 
which probe even lower radio fiuxes, also disfavor a pop- 
ulation with very low values of R below those considered 
here. 

In Figure [16] we plot our results along with two deter- 
minations of the raw obser ved R d istribution available 
in the literature , that of llvezic et al.l ((2004) and that of 
Ciras uolo et al.l ( 120061 ) along with our results. These de- 
terminations from the literature are scaled to the ratio of 
5 GHz radio to 2500 A optical luminosity assuming opti- 

Note that we have not included the density evolution which 
will shift the curves vertically but not change their shape. 
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Fig. 15. — The local distribution Gr{R) in the 5 GHz radio to 
2500 A optical luminosity ratio R for the quasars in the dataset, 
plotted as X Gj^R). The circles are from Gji>{R') as determined 
by the method of equation [S] taking account of the truncations and 
correlations in the luminosity evolutions, while the triangles result 
from forming a distribution with a naive use of the objects' raw 
ratio. The error bars represent the extremal possibilities of the 
radio-optical luminosity correlation being entirely induced (open 
circles) and entirely intrinsic (closed circles). The normalization 
is arbitrary. Also shown is the proper radio loudness distribution 
GflR, z) at redshifts z = 1 (dashed line) and z = 3 (dash-dot line), 
and z = 5 (dotted line) evolved according to the form of equation 
^ A comparison with data from the literature is shown in Figure 

cal and radio spectral indices of 0.5 and 0.6, respectively. 
In the former case (red and orange squares on the figure), 
the dataset is from an SDSS x FIRST sample which is 
fiux limited in radio and optical as in this work. The raw 
R distribution seen there is similar to the one seen here, 
although with a more pronounced dip at ~ 10. In the 
later case (blue stars on the figure), the data are from a 
small patch of sky with deep VLA observations where the 
sample of SDSS quasars has 100% completeness down to 
'S'i.4GHz = 60 /iJy, and with such radio completeness we 
would expect the observed distribution to be consistent 
with our reconstructed intrinsic distribution at low val- 
ues of which indeed it is, but not necessarily at high 
values of R, where there may be sources missing which 
have significant radio fiux but insufficient optical ffux to 
be included in the SDSS sample (see the discussion in 

8. TESTS OF ASSUMPTIONS AND ERROR ANALYSIS 

Luminosity dependent density evolution: One may be 
concerned that luminosity dependent density evolution 
(LDDE), which is not considered in the functional forms 
for the LF used here, may be necessary to represent the 
evolution of the LF. As a test of this scenario, we divide 
the data into high and low halves of de-evolved lumi- 
nosity (cutting on optical luminosity), and check the 
similarity of the computed density evolutions for the two 
sets. The density evolutions computed for both cuts are 
similar to each other, with the high half peaking in p{z) 
at z = 1.3 and the low redshift half peaking in p{z) at 
z = 1.1. Given the similarity of these distributions to 
each other and to that computed from the dataset as 
a whole, we conclude that we are justified in neglecting 
LDDE for the purposes of this analysis. 



-10 12 3 4 

Log R 

Fig. 16. — Same as FigurefTSlexcept without the evolved Gj^R, z) 
at other redshifts, and we also overpl ot the raw observed Gj^R) 
distributions obtained by llvezic etal . (2004) for 18.7 < i < 18 
(orange squares) and 15 < z < 19 (red squar es) for an SDSS x 
FIRST sample, and by Cirasuo To et al.l (|2006l ) for a deeper radio 
sample which is 100% complete to Si.4ghz = 60/iJy (blue stars). 

Optical measurement errors: There is the possibility 
that errors in the optical magnitude s could lead to a 
bias which could affect the results (jEddingtonI ll94Qf ) . 
The bias introduced by these errors would be negligi- 
ble if the source counts as a function of fiux, i.e. the 
log A/"— logs' distribution, was fiat (i.e. N{> S) oc 5^ and 
dN/dS (X S~^). Since the number density of sources in- 
creases with decreasing fiux, it is more likely that a source 
will be included than excluded. However, the magnitude 
of this effect will depend on the faint end source counts 
slope, and the shallower the slope the smaller the ef- 
fect. For constant fractional fiux measurement errors, 
an error will be introduced on the normalization of the 
source counts, and therefore the luminosity functions, 
and can be approximated by [1/2 ^heiow {"^beiow + 1)] 
(Teerikorpi 2004) where 6 is the fractional error in fiux 
and rribeiow is the faint end differential source counts 
power law slope. Typical reported SDSS DR7 quasar 
i band magnitude errors are between 0.001 and 0.002. 
However, to be conservative, we will adopt errors of 0.2 
magnitudes in i band, to account for the instrinsic RMS 
scatter due to source variability. For the faint end mag- 
nitudes of 19.1 and this error, S is around 0.16. For 
"mbeiow ^ 2 the bias on the luminosity function normal- 
ization will be around 7.6%. On the other hand, there 
will be an effect on the reconstructed slope of the lumi- 
nosity function only if the fractional measurement errors 
change with luminosity. The data show only a modest 
dependence on i band reported error with magnitude at 
magnitudes below 19.1, with magnitude errors at most 
a factor of 1.5 higher at the highest magnitude end of 
that range than for the lowest magnitudes, which corre- 
sponds to a small fractional fiux error. Because data at 
a wide range of luminosities corresponds to a given fiux, 
we consider this effect to be negligible, even if scatter due 
to variability has some fiux dependence. 

Radio incompleteness: Lastly, the selection function 
for the FIRST objects used here might not be a sharp 
Heavyside function at a peak pixel fiux density of 1 mJy, 
but rather smeared out. According to Figure 1 of JOT, the 
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selection function of FIRST for SDSS optically identified 
quasars in SDSS data release 3 is such that at an inte- 
grated flux density of 1 mJy only about 55% of sources 
are seen, and this number rises to 75% at 1.5 mJy and 
about 85% at 2 mJyQ We have considered the sample 
to be Hmited by the peak pixel flux (i.e. surface bright- 
ness limited) in the radio rather than being limited by 
the integr ated flux, in accordance with the criteria set 
forward in iBecker et al.l 0995). 

The way to test the effects on our analysis of possi- 
ble radio incompleteness is to repeat the analysis limit- 
ing the sample to a higher radio flux, where the sample 
would presumably be more complete, and determine the 
extent to which the calculated parameters change in a 
systematic way. We have done so with lower radio flux 
limit of 2 mJy (4251 objects in the canonical sample), 
as opposed to the original 1 mJy (5445 objects in the 
canonical sample). The effect, propagated through the 
analysis, is quite minor. The la uncertainties on /copt 
and /crad are extended slightly. There is no discernable 
effect on the correlation parameter a. There is a negligi- 
ble effect on the density evolutions, a small effect on the 
local luminosity functions which are reported in ^ and 
a small effect on the shape of the Gr{R) distribution 
which is sub-dominant to the effect of considering the 
radio-optical luminosity correlation to be intrinsic or in- 
duced. To the extent that faint flux radio incompleteness 
is present in the sample considered here, it does not seem 
to have a large systematic effect on the determination of 
the parameters in, and conclusions of, this analysis. 

9. DISCUSSION 

We have used a general and robust method to deter- 
mine the radio and optical luminosity evolutions simulta- 
neously for the quasars in the SDSS DR7 QSO x FIRST 
dataset, which combines 1.4 GHz radio and z-band op- 
tical data for thousands of quasars ranging in redshifts 
from 0.64 to 4.82 and over 5 orders of magnitude in radio 
loudness. 

As can be seen, the results for the SDSS DR7 QSO x 
FIRST dataset employed here are similar in bulk to the 
results for the White et al. dataset that we presented in 
QPl. In this work we show that, importantly, the major 
conclusions are not sensitive to whether the correlation 
between the radio and optical luminosities is considered 
to be induced by similar redshift evolution or intrinsic. 
In the previous work, we assumed the correlation to be 
intrinsic. Additionally, the luminosity evolution found 
for the optical only dataset matches that found for the 
combined radio-optical dataset, providing a test that the 
technique has properly handled the truncations and cor- 
relations inherent in the data. 

We have seen that the results of this analysis are nearly 
identical for the two different radio-optical matching ra- 
dius criteria considered, 1) that adopted by J07 with a 
5" radius for single matches and a 30'' radius for mul- 
tiple matches, and 2) a universal 5" matching radius 
which we feel is more appropriate. The reason for this 
preference is that 30'' radius corresponds to a physical 

The fuzzyness of the truncation boundary has a similar effect 
as the data measurement errors in the sense that it is unimportant 
for mijeiow = cind more important for larger deviations from 
this value. 



scale > 100 kpc projected^ at any redshift z > 0.2. Ex- 
tended lobes in radio galaxies are known to reach similar 
or even larger sizes. These are however objects viewed at 
large inclinations, unlike quasars obs erved, by de finition, 
at much smaller viewing ang les (Barthel'^T989^. With 
such small viewing angles, projection effects are expected 
to limit the projected sizes of the extended structures 
in radio-loud quasars. In addition, evolutionary effects 
may play a role as well, with distant quasars being typi- 
cally younger, and therefore more compact in radio, than 
nearby evolved radio galaxies. This is supported by a de- 
tailed investigation of radio morpholo g ies of SDSS (DR3) 
X FIRST quasars bv iDe Vries et al.l (|2006[ ). who noted 
that overwhelming majority of the selected sources are 
characterized by compact radio morphologies. For all 
these reasons, we believe that applying a universal 5" 
matching radius criterion in assigning FIRST counter- 
parts to the SDSS bona fide quasar sources is more ap- 
propriate, while 30" criterion for multiple matches has 
a high probability of including radio flux from objects 
unrelated to the optical quasar that they are supposedly 
linked to. 

9.1. Differential evolutions toward radio loudness with 
increasing redshift 

Quasars are seen to exhibit positive luminosity evo- 
lution in both wavebands, with stronger positive evolu- 
tion in radio than optical. The differential evolution can 
be seen in the evolution of Gr'{R') as shown in Figure 
[T5l This is in agreement with the results presented for a 
combined FIRSTxPOSS-I dataset in QPl, but in likely 
disagreement with J07 which claimed that the fraction 
of RL quasars decreases with increasing redshift. How- 
ever, Miller et al. (1990) have noted that the fractio n of 
RL quasars may increase with redshift and Donoso e t al.l 
([20091 ) reached the same conclusion by computing radio 
and o ptical LFs at different redshifts. Cirasuolo et al.l 
()2006f ) also fould that the RL fraction of quasars may 
modestly increase at high redshift. 

We attribute our apparent disgreement with the J07 
result to three causes: 1) They present one particular 
moment of the R distribution — the "radio loud frac- 
tion" — in bins, while we are presenting the distribution 
itself, so direct comparison is rather difficult, 2) they 
are calculating this fraction by including optical sources 
that have no radio detection, which is a potentially cor- 
rectable bias, and 3) they are not including the effect of 
neglecting sources that are radio bright enough to appear 
in FIRST but don't appear in the SDSS quasar sample 
that they use. 

There are potentially many AGN sources in the FRIST 
catalog that are not present in the SDSS quasar sample, 
due to their low optical fluxes. The sample used in J07 to 
calculate the radio loud fraction in bins would be missing 
this population of sources with detectable radio emission 
but lower optical fluxes than would allow an SDSS detec- 
tion or classification as a quasar source, while we have 
used an analysis which takes both flux limits properly 
into account. While for the joint dataset the present work 
has used only quasars detected in both radio and opti- 
cal and including both optical-color-selected and radio- 
match-selected sources, J07 used as their sample quasars 
detected in optical regardless of a radio detection and 
their quasars were selected based solely on optical color. 
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Stepping back, there are four basic options when con- 
structing a combined optical-radio dataset of quasars for 
an analysis of the radio loudness distribution and its evo- 
lution, given that many more confirmed quasars have an 
optical detection than a radio one: 
L quasars detected in both radio and optical 
IL quasars detected in optical regardless of a radio de- 
tection 

A. optical quasar candidates selected for followup 
based on either optical colors or presence of a radio match 

B. optical quasar candidates selected for followup 
based only on optical colors 

Under the above rubric, we have used I A for the joint 
radio-optical dataset, while JOT used IIB. 

We believe that while using option II is appropriate to 
investigate the optical luminosity function and its evolu- 
tion, such a dataset is not appropriate for investigating 
the radio loudness distribution and its evolution. This is 
because with such a dataset one will always be introduc- 
ing a bias by not including objects bright in radio but dim 
in optical, and at the same time, not being able to say if 
such a population exist at all by means of some censored 
data method. The situation could be different if one were 
able to select two quasar samples based separately on the 
radio and optical data, down to given radio and optical 
flux limits of the surveys, and then to match the two 
samples forming a master-list with both radio and opti- 
cal fluxes provided for each object included. Only then 
non-detections in either radio or optical bands could be 
considered, and the censored da ta method could be ap- 
plied (see in this c ontext, e.g., fFeigelson fc BergI 119831 : 
iFeigelso n fc NelsonI 119851 : llsobe et al.l Il986f ). However 
since an optical identification is needed to claim a quasar 
nature of a source the first place, this is not possible, so 
in using an only optically detected set one will be always 
biased against a radio-loud population. On the other 
hand, if, as in the present work, one uses option I, and 
then analyzes the data in such a way as to take the flux 
limits into account, an unbiased reconstruction of the 
intrinsic properties of the population can be achieved. 

The question of whether option A or B is most 
appropriate is less straightforward. However, quasar 
optical color may b e correlated with radio loudness 
(|Kimban et al.l 1201 If ). If this is the case, then by se- 
lecting quasars based on the optical colors solely, one 
may introduce some bias when dealing with the radio 
loudness distribution and evolution. Option lA means 
that the color criterion is not important: by including 
only those quasars which are detected in radio, and by 
allowing for radio-selected optical quasar candidates, we 
basically end up with a radio-selected (and not color- 
selected) and radio and optical flux limited sample of 
quasars. 

9.2. Lack of bi-modality in radio loudness 

Our results favor a single population of quasars with 
no bi-modality in the radio loudness parameter for the 
range of R considered here, in agreement with the con- 
clusions of QPl. Some physical implications of a single 
population in this range of R are discussed in that work, 
including that there is thus no evidence for the existence 
of two physically different population of quasars with re- 
spect to the radio/jet production efficiency. Although it 
is a somewhat different issue, there have also been studies 



of radio populations including not only quasars, but also 
Seyferts and other types of AGN in the plane of radio 
loundess versus accretion rate. Note that by definition 
all quasars accrete at high and very high rates, unlike 
radio galaxies or Seyferts. These results indicate that 
there may be two 'tracks' in the space, characterized by a 
similar non-monotonic dependence of the radio loundess 
parameter on the Eddington ratio, but differing in nor- 
malization, and corresponding to what might be consid- 
ered as an RL population hosted exclusively by elliptical 
galaxies, and an RQ population of AGN hosted by both 
types of galaxies (jSikora et al.l 120071 ). although the re- 
sult is dependent on considerations such as whether core 
or to tal radio luminosities are used (jBroderick fc Fenderl 
I2OIII ). This problem is however distinct from the one 
addressed here, deahng strictly with a well-defined pop- 
ulation of quasar-type AGN. 

In ^ we explored additional arguments against the 
presence of a bi-modality in R at lower values of R^ be- 
low those considered in this analysis. The caution to 
the above conclusion is that the sample analyzed here is 
dominated by RL objects, as we do not, for the reasons 
discussed above, consider SDSS quasars with no detected 
FIRST counterparts. Hence we cannot formally exclude 
the possibility that there may be a population of very 
radio quiet quasars which would form a separate branch 
in Figure [15] were the abscissa extended to the left, re- 
sulting in a bi-modal distribution of radio loudness for 
the entire quasar population. 

Such a numerous population of very radio-quiet 
quasars, even if absent in the FIRST database, would not 
be missed in deeper rad io surveys probing sub-mJy flux 
levels. In fact. lKimball et al. (2011) have used EVLA ob- 
servations to detect 179 SDSS DR7Q quasars in the range 
0.2 < 2; < 3 in their field down to a 6 GHz flux limit of 
20 /iJy. They claim to detect almost every quasar (97%) 
in the field in radio, and construct the 6 GHz radio lu- 
minosity function with the V/Vmax method. Their radio 
luminosity function is well fit with two populations, one 
with radio emssion dominated by AGN central engine 
emsision, and one with radio dominated by star forma- 
tion. The later becomes prominent at luminosities below 
^6GHz ^ 10^^ erg s~^ Hz~^. This corresponds to radio 
luminosities just below those probed in this analysis, as 
can be seen in Figure [lH and would imply a two pop- 
ulation population model with the more RQ population 
dominating at ~ 1 and below, although without a 'dip' 
a t ^ 10. 

i Mahonv et al.l ()2012l ) claim no bimodality in the R20 
radio loudness parameter with X-ray selecte d sample 
with z < 1 down to ~ 0.2 mJy. The iKimball et al.l 
(|2011f ) result disfavors a population with values of R be- 
low R < 0.1, finding few objects there with a volume- 
limited sample probed to low radio fluxes. Both of these 
works claim evidence for a low radio flux population 
of quasars where the radio emission is significantly en- 
hanced or even dominated by star fomation activity in 
the host galaxies. This is largely a different issue than 
the "traditional" bimodality question in quasars, i.e. as 
to whether there is a two-component shape to the radio 
loudness distribution centered around R = 10. 

The measured dN/dS distribution of extragalactic ra- 
dio sources at low fluxes, together with the limiting level 
of the extragalactic radio background (Singal et al...2010l : 
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iFixsen et al.l |2Q11[ ) and the quasar optical LF as con- 
structed here, could in principle provide further con- 
strains on the distribution of the radio-loudness param- 
eter at very low values of R. This will be addressed in a 
forthcoming analysis. 

Funding for the SDSS and SDSS-II has been pro- 
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APPENDIX 

NOTES ON THE EFFECTS OF CORRELATED RADIO AND OPTICAL LUMINOSITIES FOR THIS ANALYSIS 

As discussed in §3.2[ in the case of intrinsically correlated luminosities, the luminosity functions in different wavebands 
are not a-priori separable and one must do a coordinate rotation of the form in equation [6] to a correlation reduced 
luminosity to achieve a separable luminosity function. In that case, one must determine the two now independent 
luminosity-redshift correlation functions ^opt(^) and ^crr(^) which describe the luminosity evolutions. The full procedure 
is detailed in ^ and this result can be compared to the evolutions determined from assuming no intrinsic (i.e. only 
redshift induced) correlation between the radio and optical luminosities, to determine whether the effect is significant. 

The local LFs of uncorrelated luminosities L^p^ and L'^^^ can then be used to recover the local radio LF by a straight 
forward integration over L'^^^ and the true local optical LF as 



V^rad(i^rad) = 

As stated above this procedure can be used for the determination of the radio LF at any redshift, from which one can 
deduce that the radio luminosities also undergo luminosity evolution with 

^rad(^) = QctIz) X [^opt(^)]" (2) 

(cf equation [6|) 

Similarly we can determine the local distribution of the radio to optical luminosity ratio, R' = L'^^^/ L'^^^ = L'^^^ x 

^opt ^ ^fid "5 as 

and its evolution 

giiz) = gUz) X [goUzT-' = — (4) 

^opt 

These can then be compared with ^rad(^rad) determined from assuming no intrinsic correlation between 

the radio and optical luminosities, which is done in Figure [M] and Figure [T5l We see that the main conclusions of this 
work are not dependent on whether the observed correlation between Lrad and Lopt is intrinsic or induced. 



